Analytic solution for Brownian motion in three 
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Abstract 

Brownian motion in the three dimensional Lobachevsky space or hyperboUc space 



f~^ ' is considered in the paper written by F.I.Karpelevich, V.N.Tutubahn and M.G.Shur. 

'nJ" . A solution for radial symmetric diffusion equation in the three dimensional hyper- 

bolic space is given in that paper. However, derivation of it is not explicitly shown. 
(^ • Therefore, the diffusion equation is solved analytically with an initial condition cor- 

■^ . responding to non-zero space variable at the initial stage. 

1 

^'. 1 Introduction 

c^ ' A radial part of Brownian motion in the three dimensional Lobachevsky space or hyperbolic 

space is written by the diffusion equation, 

df D d f . ^2 df\ 

A solution of the diffusion equation is given in Ref. [1] as 



[ATcDtf sinhp ADt 

where D denotes a diffusion constant. 

We have applied Eq.(PQ) to the analyses of large transverse momentum distributions 
of secondary particles observed in high energy nucleus-nucleus (AA) collisions [2|. We 
consider p in Eq.(^ as radial rapidity, which is written with energy E, momentum p and 
mass m of observed particle by 

, ^+|P| 
p = In- 



m 
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Transverse momentum pt is defined by pt = |p| sin^, where 6 is the polar angle of the 
observed particle measured from the direction of incident colliding nuclei.^ 

However, derivation of solution (^ is not explicitly shown in Ref . jT] . In addition, up 
to the present, we cannot find the article where the analytical proof is made.^ Therefore, 
we would like to solve the diffusion equation analytically. 

2 Diffusion equation in three dimensional hyperbolic 
space 

The radial symmetric diffusion equation in the three dimensional hyperbolic space H^ is 
given by. 

The initial condition for Eq.(j2I) is taken as, 

/(,,«.( = 0) = ^<^. (3) 

47rsmh p 

At first, f{p,po,t) is put as 

f{p,Po,t) = -^e-'''g{p,t). (4) 

smhp 

From Eqs.Q and (JD), equation for g(p,t) is given by 



dt dp 



- D^2 = 0' (5) 



which is the diffusion equation in one dimension. The initial condition for Eq.(j3)) is given 
from Eqs.Q and (jlj as. 



g{p^t = 0) = - — — • (6) 

Att smh p 



The Laplace transform. 



c» 



G{p, s) = C[g{p, t),s]= g{p, t)e-''dt, (7) 

Jo 



^In high energy nucleus-nucleus collisions, colliding energy y^ grows up to -^/i — 200 GeV, and sec- 
ondary particles with transverse momentum pT more than 10 GeV/c are observed. The observed pr 
distributions have long tail compared with exponential distribution in px- As is well known, the funda- 
mental solution of stochastic process is Gaussian if variables in the euclidian space are used. Therefore, 
as long as we consider stochastic equations in the transverse momentum space, it is very hard to describe 
observed pt distributions. This fact suggests that a relativistic approach to the stochastic process would 
be needed. 

^In Ref. |3] , the solution is obtained from a Markovian property of Brownian motion. 



is applied to Eq.©; 



kG{p,s) = —-g{p,t = 0), 



dp^ '^' ' D 



K = ^s/D. (8) 

In order to solve Eq.ijSJ, we put 

G{p,s) = C{p)e-'^^. (9) 

From Eqs.dHl) and Q, we obtain 

Equation (fTUI) is rewritten as, 

After integrating over p, we obtain, 
dC{p) 



u Jo 



dp D Jo 

where -B(O) = dC{0)/dp. Then, C{p) is given in the following form, 

C(P) = C'(O) + ^ (e'-" - l) - 7^ re'^-^'dp' f g{i,t = 0)e~-^d^, 
2k ^ ^ D Jo Jo 

= c(0) + ^ (e'^'' - l) - ^ [ 3{i. t = 0) {e'-P-^^ - e<) d^ (11) 

If p > po, after substituting the initial condition (jH)) to Eq. (fTT|l . we have an expression for 

C(P); 

C(rt = C(0, + - (b(0) - j^^j-j^j + - (j^^^-j^ - B(0)] . 

Then, ^(p, s), the Laplace transform of g{p,t) is written as 
Gip,s) = Cip)e-^^ 

ftp / p-ftPO \ p-Kp / pKpo \ 

= C(0)e-'''' + — 5(0)-— ^ ^ + ^ ( J^ , 5(0) . 

^ ^ 2fi; V ^ ^ 47rDsinhpo/ 2k ^vrDsinhpo ^ 7 

In order to get inverse Laplace transform of G{p, s), equation, 

5(0) 
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should be satisfied. Then, we have 



G{p, s) = C(0)e-- + 1 ^inhj.po ^-.p^ ^^2) 

AtiDk smhpo 



Substituting k = Js/D to Eq. (jl2j) . we obtain, 



1 1 ^ / s\"' pI"'^^ 



The inverse Laplace transform of G{p, s) is defined by, 

]^ pc+ioo 

g{p,t) = C-^[G{p,s);t] = 7— G{p, s) exp[st\ds, 

Zm Jc-ioo 

where c is a constant. The following formulae for the inverse Laplace transform^. 



£-^[exp[-pv^//D];t] = ^ ^ exp 



P „.„r P 
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are useful. In the above equation, Hn{x) denotes the Hermite polynomial of degree n. 
Then, we find 



{ATiDty 4£)t 47rVvr-Dtsinhpo 

00 -1 / \ 2m+l 2 



3 Solution for radial symmetric diffusion equation 

From the generating function [64 of the Hermite polynomial. 



2m+l y.2m+l 

e"*' sinh 2a;t = JZ 77^ ^T7^2m+i(a;), 

^0 (2m +1)! 



and Eq.(fT^. we find 



/(p, Po, t) = -Ve-^*(7(p, t) = 47rDC(0)/i(p, t) + Mp, po, t), 
smhp 

/i(P,t) = ,, ^ , e-^^-4-exp[ ^'1 



(47rDt)3 sinhp 4Dt ' 
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Now, we extend the definition region of f{p,po,t) from p > po to p > 0. Function 
f2{p, Poyt) satisfies 



/ dp I dO I d(f) f 2{ p, po,t) smh."^ p sin 6' = 1. 
Jo Jo Jo 



dp de 

10 Jo Jo 

Function f{p,po,t) is tlie probability density, and satisfies tlie normalization condition, 

/•oo rn t'2'K 

/ dp dO d(f) f {p, po,t) smh"^ p sin6' = 1, 
Jo Jo Jo 

from which constant C(0) is determined as 

C(0) = 0. 
Then, the solution of Eq.Q, which satisfies the initial condition Q, is given by 

f(nn f)- ^ c-^' si^h(lz^) ,,.^r P' + Po i ^ S^ 

j[P^Po,t - - — , e -T—r r-r^exp „ ■ [i-O 

27rV47rDt smhposmhp ADt 

4 Summary and discussions 

In the present paper, we have investigated the radial symmetric diffusion equation Q in 
the three dimensional hyperbolic space H^ with the initial condition Q, which corresponds 
to the non-zero space value at the initial stage. The solution is give by Eq.(|15|). 

It should be noted that Eq.(fT3j) reduces to Eq.(IT]) in the limit of po ~^ 0. Namely, 

f{p,Po,t) satisfies, 

lim f(p,Po,i) = f{p,t). 

Po^O 

It will be necessary to consider a physical meaning of po iii Eq. ()15|) . In high energy 
AA collisions, some nucleons inside the nucleus would suffer inelastic collisions more than 
once. In such a case, colliding nucleons have some momentum distribution in the center of 
mass system. Therefore, secondary particle system produced in a nucleon-nucleon collision 
would have non-zero velocity even at the initial stage in the center of mass system. Such 
a situation would correspond to the initial condition Q. Therefore, a sort of velocity 
distribution of secondary particle systems at the initial stage is included in Eq. (fT3j) . 

The radial symmetric diffusion equation in the Euclidian space corresponds to the 
Bessel process jHl El- It is generalized to the Bessel process X^"''^^ in the wide sense with 
index (a, c) in Ref. |3]. The process X^°''^^ is the stochastic process in the «-dimensional 
hyperbolic space Ha- The transition density of X^"''^^ is denoted by p^^'^^t, x, y) in Ref. jS], 
and derived from a probabilistic approach different from ours. In fact, the solution (fT3j) 
for the diffusion equation (J2]) with the initial condition satisfies. 

However, relation between p and po does not appear explicitly in Ref. [3]. 
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